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Abstract. Massive gravity theories have been developed as viable IR modihcations of gravity 
motivated by dark energy and the problem of the cosmological constant. On the other hand, 
modihed gravity and modihed dark matter theories were developed with the aim of solving 
the problems of standard cold dark matter at galactic scales. Here we propose to adapt the 
framework of ghost-free massive bigravity theories to reformulate the problem of dark matter 
at galactic scales. We investigate a promising alternative to dark matter called dipolar dark 
matter (DDM) in which two different species of dark matter are separately coupled to the 
two metrics of bigravity and are linked together by an internal vector held. We show that 
this model successfully reproduces the phenomenology of dark matter at galactic scales {i.e. 
MOND) as a result of a mechanism of gravitational polarisation. The model is safe in the 
gravitational sector, but because of the particular couplings of the matter helds and vector 
held to the metrics, a ghost in the decoupling limit is present in the dark matter sector. 
However, it might be possible to push the mass of the ghost beyond the strong coupling scale 
by an appropriate choice of the parameters of the model. Crucial questions to address in 
future work are the exact mass of the ghost, and the cosmological implications of the model. 
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1 Introduction 

In the last century, cosmology has progressively developed from a philosophical to an empirical 
scientihc discipline, witnessing high precision cosmological observations, which culminated 
with the standard model of cosmology, the A-CDM model [1]. The standard model is based 
on General Relativity (GR) and is in great agreement with the abundance of the light elements 
in the big bang nucleosynthesis, the anisotropies of the cosmic microwave background (GMB), 
the baryon acoustic oscillations, lensing and the observed large scale structures. It notably 
relies on the presence of dark energy in form of a cosmological constant A, giving rise to the 
accelerated expansion of the universe. 

1.1 Massive gravity context 

In the prevailing view, the cosmological constant corresponds to the constant vacuum energy 
density which receives large quantum corrections. Unfortunately, so far there is no successful 
mechanism to explain the observed unnatural tiny value of the cosmological constant, see 
e.g. [2, 3]. This difficulty has sparked a whole industry studying modihcations of gravity in 
the infra-red (IR) invoking new dynamical degrees of freedom. 

One tempting route is massive gravity, motivated by the possibility that the graviton 
has a mass. It was a challenge over forty years to construct a covariant non-linear theory for 
massive gravity. The foregathered remarkable amount of effort hnally rose to the challenge [4- 
11], to construct the potential interactions in a way that got rid of the Boulware-Deser (BD) 
ghost [12]. This theory has been further extended to more general models by adding additional 
degrees of freedom. A very important extension is massive bigravity theory [13]. 

Once the interactions in the gravitational sector were guaranteed to be ghost-free, a 
natural follow up question was how to couple this theory to the matter helds without spoiling 
the ghost-freedom. First attempts were already discussed in the original paper of bigravity [13]. 
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If one couples the matter fields to both metrics simultaneously, this reintroduces the BD 
ghost [14, 15]. Furthermore, the one-loop quantum corrections detune the special potential 
interactions at an unacceptable low scale and hence this way of coupling is not a consistent 
one [15]. The safer way is to couple the matter helds to just one metric. In this way, the 
quantum corrections give rise to contributions in form of cosmological constants. One could 
also try to couple the matter held to the massless mode, which is unfortunately also not ghost 
free [16]. Another possible way of coupling the matter helds to both metrics is through a 
composite ehective metric constructed out of both metrics [15, 17-19], which is unique in the 
sense that it is the only non-minimal matter coupling that maintains ghost-freedom in the 
decoupling limit [20-22]. Furthermore, the quantum corrections are guaranteed to maintain 
the nice potential structure. Other important consequence of this new matter coupling is the 
fact that it helps to evade the no-go result [23] for the hat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background. A detailed perturbed ADM analysis revealed the existence of a 
BD ghost originating from an operator involving spatial derivatives [15, 17]. Therefore, the 
ghost will probably reappear for highly anisotropic solutions. Since the ghost remains absent 
up to the strong coupling scale, the matter coupling can be considered in an ehective held 
theory sense at the very least till the strong coupling scale. The precise cut-oh of the theory 
or mass of the ghost has still to be established. 

The absence of the BD ghost is not only important at the classical level, but also at the 
quantum level. For the classical ghost-freedom the relative tuning of the potential interactions 
is the key point. Therefore one has to ensure that the quantum corrections do not detune 
the potential interactions. Concerning the decoupling limit, it is easy to show that the theory 
receives no quantum corrections via the non-renormalization theorem due to the antisymmetric 
structure of the interactions [24] (in fact the same antisymmetric structure of the Galileon 
interactions protect them from quantum corrections [25-27]). Beyond the decoupling limit, 
the quantum corrections of the matter loops maintain the potential interactions provided that 
the above criteria are fulhlled. Concerning the graviton loops, they do destroy the relative 
tuning of the potential interactions. Nevertheless, this detuning is harmless since the mass of 
the corresponding BD ghost is never below the cut-off scale of the theory [28]. The bimetric 
version of the theory shares the same property [29]. 

Massive gravity/bigravity theory has a rich phenomenology. The decoupling limit admits 
stable self-accelerating solutions [30], where the helicity-O degree of freedom of the massive 
graviton plays the role of a condensate whose energy density sources self-acceleration. Unfor¬ 
tunately these solutions suffer from strong coupling issues due to the vanishing kinetic term of 
the vector modes [30, 31]. With the original massive gravity and the restriction of flat reference 
metric, one faces the no-go result, namely that there are no flat FLRW solution [23]. One can 
construct self-accelerating open FLRW solutions [32], which however have three instantaneous 
modes [33] and suffer from a nonlinear ghost instability [34]. The attempt to promote the ref¬ 
erence metric to de Sitter [35, 36] also failed due to the presence of the Higuchi ghost [35]. In 
order to avoid these difficulties, one either gives up on the FLRW symmetries [23], or invokes 
new additional degrees of freedom [13, 37, 38]. 

Thanks to the freedom gained in the inclusion of the second kinetic term in the bimetric 
extension [13], there now exists many elaborate works concerning the cosmology of the bigrav- 
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ity theory, see e.g. [39-41]. In the case of minimally coupled matter helds and small graviton 
mass, the theory admits several interesting branches of solutions. Unfortunately, among them 
the self-accelerating branch is unstable due to the presence of three instantaneous modes, 
and a second branch of solutions admits an early time gradient instability [42]. Nevertheless, 
there exists attempts to overcome the gradient instability either by viable though hnely tuned 
solutions in the case of a strongly interacting bimetric theory with m ^ Ho [43, 44], or by 
demanding that Mg 3> Mf as was proposed in [45]. See also other recent works concerning 
the phenomenology of bimetric gravity [46-50]. 

1.2 Dark matter context 

On a quite different but equally fascinating topic is the phenomenology of dark matter at 
galactic scales. The evidence for dark matter is through the measurement of the rotation 
curves of spiral galaxies which turn out to be approximately flat, contrary to the Newtonian 
prediction based on ordinary baryonic matter [51, 52]. The standard explanation is that the 
disk of galaxies is embedded into the quasi-spherical potential of a huge halo of dark matter, 
and that this dark matter is the same as the cold dark matter (CDM) which is evidenced 
at large cosmological scales notably with the £t of A-CDM to the CMB anisotropies [53, 
54]. Unfortunately this explanation faces severe challenges when compared to observations 
at galactic scales [55, 56]. There are predictions of the A-CDM model that are not observed, 
like the phase-space correlation of galaxy satellites and the generic formation of dark matter 
cusps in the central regions of galaxies. Even worse, there are observations which are not 
predicted by A-CDM, such as the tight correlation between the mass discrepancy (luminous 
vs. dynamical mass) which measures the presence of dark matter and the involved scale of 
acceleration, the famous baryonic Tully-Fisher (BTF) relation for spiral galaxies [57-59], and 
its equivalent for elliptical galaxies, the Faber-Jackson relation [60]. 

Instead of additional, non-visible mass, Milgrom [61-63] proposed an amendment to 
the Newtonian laws of motion in order to account for the phenomenology of dark matter 
in galaxies, dubbed MOND for MOdihed Newtonian Dynamics. According to the MOND 
hypothesis the change has a relevant influence on the motion only for very small accelerations, 
as they occur at astronomical scales, below the critical value Oq — 1.2 x lO^^'^m/s^ measured 
for instance from the BTF relation.^ A more elaborate version of MOND is the Bekenstein- 
Milgrom [64] modification of the Poisson equation of Newtonian gravity. All the challenges 
met by A-CDM at galactic scales are then solved — sometimes with incredible success — 
by the MOND formula [55, 56]. Unfortunately, MOND faces difficulties in explaining the 
dark matter distribution at the larger scales of galaxy clusters [65-69]. It is also severely 
constrained by observations in the solar system [70, 71]. 

Reconciling A-CDM at cosmological scales and MOND at galactic scales into a single 
relativistic theory is a great challenge. There has been extensive works on this. One approach 
is to modify gravity by invoking new gravitational degrees of freedom without the presence of 
dark matter [72-81]. Primary examples are the tensor-vector-scalar (TeVeS) theory [73, 74] 
and generalized Einstein-iEther theories [75, 76]. Another approach is a new form of dark 

striking observation is the numerical coincidence that op ~ -x/A [55, 56]. 
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matter a la MOND, called dipolar dark matter (DDM). It is based on the dielectric analogy 
of MOND [82] — a remarkable property of MOND with possible far-reaching implications. In¬ 
deed the MOND equation represents exactly the gravitational analogue (in the non-relativistic 
limit) of the Gauss equation of electrostatics modihed by polarisation effects in non-linear di¬ 
electric media. Some early realizations of this analogy were proposed in [83-85] and shown 
to also reproduce the cosmological model A-CDM. The best way to interpret this property 
is by a mechanism of gravitational polarisation, involving two different species of particles 
coupled to different gravitational potentials. This was the motivation for pushing forward a 
more sophisticated model in the context of a bimetric extension of GR [86, 87]. In this model 
the two species of dark matter particles that couple to the two metrics are linked through an 
internal vector held. The model [87] is able to recover MOND at galactic scales and the cos¬ 
mological model A-GDM at large scales. Furthermore the post-Newtonian parameters (PPN) 
in the solar system are the same as those of GR. Unfortunately, even if this model delivers a 
very interesting phenomenology, it is plagued by ghost instabilities in the gravitational sector 
already at the linear order of perturbations and within the decoupling limit [88]. A more 
promising model was then proposed in [88], with the hope to cure these instabilities in the 
gravitational sector. 

In the present work we will study in detail the theoretical and phenomenological con¬ 
sistency of this recently proposed model for dipolar dark matter in the context of bimetric 
gravity [88]. We will hrst work out the covariant held equations in Sec. 2 and analyze the 
linear held equations in Sec. 3. As next, we will investigate in Sec. 4 the mechanism of gravita¬ 
tional polarisation in the non relativistic approximation, and show how it successfully recovers 
the MOND phenomenology on galactic scales. Finally we will pay attention in Sec. 5 to the 
matter sector and investigate the interactions in the decoupling limit. The outcome of Sec. 5 
is that a ghost instability is still present in the (dark) matter sector of the model. The paper 
ends with some short conclusions in Sec. 6. 

2 The covariant theory 

A new relativistic model for dipolar dark matter has been recently proposed in [88]. Basically 
this model is dehned by combining the specihc dark matter sector of a previous model [87] 
with gravity in the form of ghost-free bimetric extensions of GR [13]. The dark matter in this 
model consists of two types of particles respectively coupled to the two dynamical metrics 
and of bigravity. In addition, these two sectors are linked together via an additional 
internal held in the form of an Abelian U{1) vector held This vector held is coupled to 
the ehective composite metric Qeg of bigravity which is built out of the two metrics g and 
/, and which has been shown to be allowed in the matter Lagrangian in the ehective held 
theory sense at least up to the strong coupling scale [15, 17]. The total matter-plus-gravity 
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Lagrangian reads^ 


^_/M2 

^ V 9 ( 2 ^9 Pbar Pg 


~ Pf 


+ V~9eS 


+ A,{f^ - j^) + XM^^W{X) 


( 2 . 1 ) 


Here Rg and Rf are the Ricci scalars of the two metrics, and the scalar energy densities of the 
ordinary matter modelled simply by pressureless baryons,^ and the two species of pressureless 
dark matter particles are denoted by pbar, Pg and pj respectively. In addition j^, stand for 
the mass currents of the dark matter, dehned by 

= and = ( 2 , 2 ) 

V ^?efr V 9eS 

where = VgPgU^ and Jj = VfPfU^^ are the corresponding conserved dark matter currents 
associated with the respective metrics g and /, thus satisfying '^gJg = 0 and '^gJf = 0. 
Here Pg and pf denote the scalar densities and are the four velocities normalized to 

PguUgUg = — 1 and fgiyU^u’j = — 1. We also introduced two constants Vg and Vf specifying the 
ratios between the charge and the inertial mass of the dark matter particles. 

The metric independent matter degrees of freedom are the coordinate densities p* = 
yJ—gPgU^g and p*f = y^—fpfU^f, and the coordinate velocities = Ug/u^ and Vf = u^f/vPf. 
The associated metric independent currents are J*^ = rgp*v'^ and = rfp*fVf. They are 
conserved in the ordinary sense, d^J*^ = 0 and = 0. They relate to the classical 

currents J^, Jj or j^, by 




Notice that the coupling term y/—g^Ag{i^ ~ action (2.1) is actually independent 

of any metric, neither g nor / nor ges- We shall study in detail in Sec. 5 the implications of 
the term — jf) for the decoupling limit of the theory. 

The vector held Ag, is generated by the dark matter currents and plays the role of a 
“graviphoton”. The presence of this internal held is necessary to stabilize the dipolar medium 
and will yield the wanted mechanism of gravitational polarisation, as we shall see in Sec. 4. 
It has a non-canonical kinetic term W(<T), where W is a function which is left unspecihed at 
this stage, and 


T’ 


-pgu -p 

IIU 

4A 


(2.3) 


Here A is a constant and the held strength is constructed with the ehective composite metric 
given by (2.4) below, i.e. = g^Pla^pa with = dgAy - d^Ag not depending on 
the metric as usual. 


^The metric signature convention is (—, +, +, +). We adopt geometrical units with the gravitational con¬ 
stant, the Planck constant and the speed of light being unity, G = h = c = 1, unless specified otherwise. 

^We have in mind that the baryons actually represent the full standard model of particle physics. 


5 















The model (2.1) is defined by several constant parameters, the conpling constants Mg, 
Mf and the mass of the graviton m, the charge to mass ratios Vg and nj of dark matter, 
the constant A associated with the vector field, and the arbitrary constants a and {3 entering 
the effective metric (2.4). Parts of these constants, as well as the precise form of the fnnction 
>V(T’), will be determined in Sec. 4 when we demand that the model reprodnces the reqnired 
physics of dark matter at galactic scales. In particnlar A will be related to the MOND 
acceleration scale Oq, and the model will finally depend on oq and the mass of the graviton 
m. Additionally, there will be still some remaining freedom in the parameters a and (3, and 
in the conpling constants Mg and Mf. 

In the model (2.1) the ghost-free potential interactions between the two metrics g and / 
take the particnlar form of the sqnare root of the determinant of the effective metric [15, 17, 19] 

C = “"Sc- + 2a/3 e" + , (2.4) 

where a and f3 are arbitrary constants, and we have defined 

C = = f„Yf, (2.5) 

where the sqnare root matrix is defined hy X = together with its inverse Y = \l f~^g. 

Interchanging the two metrics g and / does not change the form of the metric (2.4)-(2.5) except 
for a redefinition of the parameters a and (3. Notice that can be proved to be antomatically 
symmetric [16, 89], and in [88] it was shown, that it corresponds to the composite metric 
considered in the previous model [87]. The square root of the determinant of gj^ is given in 


either forms respectively associated with the g or f metrics, by 

v'-fi'eff = det(a!l i3X) = det(/91 -1- aY) . (2.6) 

More explicitly, introducing the elementary symmetric polynomials e„(X) and en{Y) of the 
square root matrices X or Y,'^ 

eo{X) = l, (2.7a) 

efX) = p] . (2.7b) 

e 2 (.Y) = l(p]"-p=]), (2.7c) 

e,{X) = -3[A'] [A'2] +2[A=>]) , (2.7d) 

e4(A) = 4 ([A] ‘ - 6 [A]" [A=] + 3 [A^] V 8 [A] [A’] - 6 [A‘]), (2.7e) 

we can write, still in symmetric guise, 

4 4 

o‘‘"”/3"e>.(A) = a''P‘~"'en{Y ). (2.8) 

n=0 n=0 


"‘As usual we denote the traces of rank-2 tensors as = [Al], = [Al^], and so on. Recall in 

particular that y/~9 ^n{X) = \/~f S 4 -n(X), with e^^X) = det(N) and e 4 (F) = det(F). 
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In this form we see that (2.6) corresponds to the right form of the acceptable potential inter¬ 
actions between the metrics g and /. 

We can hrst vary the Lagrangian (2.1) with respect to the two metrics g and /, which 
yields the following two covariant Einstein held equations [90] 

3 

A/| yy G;)" - V/qy' 

n=0 

=yy { k 1 ;+ ry")+(“b" ry.'+ utz ) , ( 2 . 9 a) 

3 

n=0 

= xyr;>' + (/32fyry.'’ + ary.). (2.9b) 

where and are the Einstein tensors for the q and f metrics. The tensors U'fK and 

9 1 J {n) 

are dehned from the following sums of matrices [8]^ 

n 

^n) = 5^(-)"e„-p(X)X^ 

p=Q 
n 

VH = J2^-ren-r{Y)Y\ 

p=0 

by raising or lowering indices with their respective metrics, thus and = 

By the same property which makes (2.5) to be symmetric, one can show that 
and are indeed automatically symmetric. 

In the right sides of (2.9) the stress-energy tensors and are dehned with respect 
to the metric g, is dehned with respect to / and with respect to ges. Thus, for 
pressureless baryonic and dark matter huids we have = PbarnOa^^nbar’ '^g'" ~ Pg'^g'^'g and 

= pfU^u^j- in terms of the corresponding scalar densities and normalized four velocities, 
while represents in fact the stress-energy tensor of the vector held A^, i.e. 


(2.10a) 

(2.10b) 


TgZ = 




fiiy 

eff 


( 2 . 11 ) 


where kV;^ is the derivative of TV with respect to its argument A dehned by (2.3). 

The equation of motion for the baryons follows a geodesic law a])®''' = = 0, 

whereas the equations of motion for the dark matter particles are non-geodesic. 


Op =, (2.12a) 

al = (2.12b) 

®Here we adopt a slight change of notation with respect to [8]. Our notation is related to the one of [8] by 
t/(„) = {-rY^n){X) and E(„) = (-)-y(„)(r). 


7 









where the accelerations and af^ = u'j'Vluj^ are dehned with their respective 

metrics. On the other hand, the eqnation for the vector held yields 


ySeff 








(2.13) 


with denoting the covariant derivative associated with g^s- This eqnation is obviously 
compatible with the conservation of the currents, since by (2.2) we have also = 0 and 

= 0- The equation (2.13) can be also written as 

= (2.14) 


Combining the equations of motion (2.12) with (2.14) we obtain the conservation law 


x/=9 + ^/=7vpl„ + = o. ( 2 . is) 


Alternatively, such global conservation law can be obtained from the scalarity of the total 
matter action under general diffeomorphisms. 


3 Linear perturbations 

In this section we will be interested in computing the linear perturbations of our model, which 
will be extensively used in the following section for the study of the polarisation mechanism. 
Following the ideas of [87] we will assume that the two fluids of dark matter particles are 
slightly displaced from the equilibrium conhguration by displacement vectors and y^. This 
makes the dark matter medium to act as an analogue of a relativistic plasma in electro¬ 
magnetism. The dark matter currents will be slightly displaced from the equilibrium current 
Jo = PqUq as well, with satisfying = 0. Furthermore, we will perturb the two metrics 

around a general background in the following way® 

Pfj.i' T hfiu) 1 (3.1a) 

ffiu = + ifiuY , (3.1b) 

with main interest in background solutions where In the latter case the potential 

interactions take the form 

4 

X](a + PY~'^en{k ), (3.2) 

n=0 

where = ahfj_u -|- and traces are dehned with the common background metric g = f ■ 
Working at hrst order in the displacement of the dark matter particles and assuming that 
their gradients are of the same order as the metric perturbations, the dark matter currents 
can then be expanded as [87] 

j; = jS + 2Vr(jtj-?' <) + 0(2), 
f, = jS + 2V"(jl‘'±jl sj) + 0(2), 


®By which we mean, for instance, = {g^^p -I- hpp)gP'^{g^y + h„y). 


(3.3a) 

(3.3b) 








with the projector operator perpendicular to the four velocity of the equilibrium Mq given as 
Oes +'*^o'^ 0 ) ^he remainder 0{2) indicating that the expansion is at hrst order. The 
desired plasma-like solution is achieved after plugging our ansatz (3.3) into the equations of 
the vector held (2.13), which results in 

j'iT’ + 0(2), (3.4) 

where we have dehned the projected relative displacement as =T0 {i/g — y'j). Note that 
is necessarily a space-like vector. In Sec. 4 we shall see that the spatial components of this 
vector, which can be called a dipole moment, dehne in the non-relativistic limit the polarisation 
held of the dark matter medium as P* = where Pq is the coordinate density associated 
with pq. But for the moment we only need to notice that is a hrst order quantity, therefore 
the held strength is itself a hrst order quantity, and hence the stress-energy tensor (2.11) 
with respect to peff is already of second order,' 


TZ, = 0(2). (3.5) 

For this reason, in the case of our desired plasma-like solution the contributions of the self¬ 
interactions of the internal vector held will be at least third order in perturbation at the level 
of the action. For more detail see the comprehensive derivations in [87]. 

Let us hrst expand the Lagrangian (2.1) to hrst order in perturbation around a common 
background g = f ■ We assume that there is no matter in the background, so the matter 
interactions do not contribute at that order. Using (3.2) we obtain (modulo a total divergence) 




Z^j\-(Mpr + M]fr)G% + m=Afi(a + /?)*[«:] 




where is the Einstein tensor for the background metric p, and Cg is the background value 
of the Lagrangian. In order for the background in our interest to be a solution of the 
theory, we have to impose that linear perturbations vanish. We hnd 

= w?^(5{a + I3fgg^ . (3.7) 

This is only possible if 

With this choice we guarantee that we are expanding around the background that 

is a solution to the background equations of motion. 

As next we shall compute the action to second order in perturbations. Let us emphasize 
again that thanks to the plasma-like ansatz (3.3) the self-interaction of the internal vector held 

^We are anticipating the form of W given by (4.23) and which implies W = 0{2) and Wx = 1-1- 0{1). 
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does not contribute to this order, see (3.5). Our quadratic Lagrangian above the background 
simply becomes 


C — Cg + 


+ v(rc + ry) + ^rf 


-m^M^g(a + l3 Y{[k] - [fc^]) 


+ 0(3), 


(3.9) 


where is the Lichnerowicz operator on the background g given in the general case by 

- S,Ji) + K^lh - 

- 2C5,„ft'>' - Irs (ft"" - , (3.10) 

where h = [h] = with denoting the Weyl curvature of the back¬ 

ground metric. Of course this expression is to be simplihed using the background equa¬ 
tions (3.7). Also bear in mind that aMj = (3Mg for having a common background g = f. 
The held equations obtained by varying with respect to and yield 


- + TZ + TZ + a{a + {g^'^k - = 0(2), (3.11a) 

- = 0(2), (3.11b) 


where we recall that = ah^'^ + and k = [k] = ggyk^'^. Of course those equations can 
also be recovered directly from the general Einstein held equations (2.9). 

In the present paper we shall mostly make use of that combination of the two Einstein 
held equations (3.11) which corresponds to the propagation of a massless spin-2 held. We 
subtract the two equations from each other so as to cancel the mass term, resulting in 


fMf 


M? \ 1 1 

j + - [TZ, + TZ) - = 0(2). (3.12) 


However we shall also use the contribution of the mass term in the linearised Bianchi identities 
associated with (3.11). Thus, last but not least we can act on (3.11) with V^, yielding 


= 0 ( 2 ). 


(3.13) 


In deriving this relation we use the result that, for instance, V^TZ — PgZ ~ ^(2) which 
comes from the equations of motion (2.12) and the fact that there is no matter in the back¬ 
ground. We shall see in the next section how important is this relation for the present model 
to recover the looked-for phenomenology of dark matter at galactic scales. 


4 Polarisation mechanism and MOND 

We next consider the Newtonian or non-relativistic (NR) limit of our model. In the previous 
section we derived the held equations at linear order around a vacuum background metric g, 
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which necessarily obeys the equations (3.7). Thus, the background is a de Sitter solution, with 
cosmological constant given by the graviton’s mass, A ~ When computing the Newtonian 
limit applied to describe the physics of the local universe, e.g. the solar system or a galaxy at 
low redshift, we shall neglect the cosmological constant and shall approximate the de Sitter 
metric ^ by a flat Minkowski background rj. 

The best way to implement (and justify) this approximation is to perform an expansion 
on small scales, say r —0. When we solve for the gravitational field in the solar system or 
a galaxy embedded in the de Sitter background, we get terms like 1 — 2Mr~^ — Ar^/3 (de 
Sitter-Schwarzschild solution), where M is the mass of the Sun or the galaxy. On small scales 
the term 2Mr~^ will dominate over Ar^/3, so we can ignore the influence of the cosmological 
constant if the size of the system is -C (6M/A)^/^. In the present case we shall require that the 
size of our system includes the very weak field region far from the system where the MOND 
formula applies. In that case the relevant scale is tq = M/a^, where oq — 1.2 x 10“^° m/s^ is 

the MOND acceleration. So our approximation will make sense provided that tq (6M/A)^/^. 
For a galaxy with mass M ~ 10^^ Mq the MOND transition radius is tq ~ 10 kpc. With 
the value of the cosmological constant A (3Gpc)-2 we hnd (6M/A)V3 700 kpc, so we 

are legitimate to neglect the influence of the cosmological constant. For the solar system, 
ro ~ 0.04 pc while (OMq/A)^/^ 150 pc and the approximation is even better. 

Actually we shall make the approximation g — r] only on the particular massless spin-2 
combination of the two metrics given by (3.12), namely 

- - -^er) = (r^ + r^-) + -ry + 0(2), ( 4 .i) 

where the Lichnerowicz operator (3.10) is now the flat one, and also in the constraint equa¬ 
tion (3.13). We shall not need to consider the other, independent (massive) combination of 
the two metrics. In the NR limit when c —?> 00 we parametrize the two metric perturbations 
and by single Newtonian potentials Ug and Uf such that 

hoi — O (^^^ ’ 

y- = + o (; 

We shall especially pay attention to the potential Ug since it represents the ordinary Newtonian 
potential felt by ordinary baryonic matter. Similarly we parametrize the internal vector held 
Ag by a single Coulomb scalar potential 0 such that 


4 = 0 A) , (4.2b) 

y) > + ° (?) ■ 


A = I + O (4) . (4.3a) 

A = o (4) . (4.3b) 


11 



It is now straightforward to show that (4.1) reduces in the NR limit to a single scalar equation 
for a combination of the Newtonian potentials Ug and 


A 



/3 



(pLr + P*g) + -^P*f , 


(4.4) 


where A = is the ordinary Laplacian, and Pg P} ^^e ordinary Newtonian 
(coordinate) densities, satisfying usual continuity equations such as dtp* + V • {PgVg) = 0. 
Notably, all non-linear corrections 0{2) in (4.1) are negligible. In principle, we should a 
priori allow some parametrized post-Newtonian (PPN) coefficients in the spatial components 
ij of the metrics (4.2), say and 7/. However, when plugging (4.2) into (4.1) we determine 
from the ij components of (4.1) that in fact 7^ = 1 and 7/ = 1. Similarly we hnd that the 
equations (2.13) governing the internal vector held reduce to a single Coulomb type equation. 


V 




Jp{^9p*g-^fP*f) 


(4.5) 


Notice that the constants a and fj in the effective metric peff to which is coupled the internal 
held intervene in the NR limit only in the expression 


2\{a + (3Y ■ 


(4.6) 


Next, the equations of the baryons which are geodesic in the metric g reduce to 


dUbar 

dt 


VU, 


g ’ 


while the equations of motion (2.12) of the dark matter particles become 


(4.7) 


+ (4.8a) 

^ = Vt/, - T,Vt> . (4.8b) 

At this stage the important point is to recall that the two sectors associated with the 
metrics g and / do not evolve independently but are linked together by the constraint equa¬ 
tion (3.13). We now show that this constraint provides a mechanism of gravitational polar¬ 
isation in the NR limit and yields the MOND equation for the potential Ug felt by baryons. 
In hat space-time this constraint is 


54F" - = 0(2), (4.9) 

where and k = p'^‘^kp„. Plugging (4.2) into (4.9) we readily obtain a 

constraint on the gravitational forces felt by the DM particles, namely 

V{aUg + pUf) =d. (4.10) 

®We shall adopt the usual boldface notation for ordinary three-dimensional Euclidean vectors. Also, from 
now on we no longer write the neglected remainders D(l/c^). 
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Note that the constraint (4.10) comes from a combination between the 00 and ij components 
of the metrics (4.2). Using (4.10) we express the eqnations of motion (4.8) solely in terms of 
the potential Ug ruling the baryons, 


dt 

dt 


VUg + rgVcj), 

a 


(4.11a) 

(4.11b) 


As we see, the gravitational to inertial mass ratio mg/m\ of the / particles when measured 
with respect to the g metric is m^/mi = —a/(3. 

We look for explicit solutions of (4.11) in the form of plasma-like oscillations around 
some equilibrium solution. We immediately see from (4.11) the possibility of an equilibrium 
(he. for which dvg/dt = dvf/dt = 0) when we have the following relation between constants, 


^ = !1 
/3 r, ' 


(4.12) 


The equilibrium holds when -|- Vg'Vcj) = 0, he. when the Coulomb force annihilates the 
gravitational force. To describe in the proper way the two DM fluids near or at equilibrium 
we use the NR limits of the relations (3.3), with an appropriate choice of the equilibrium 
conhguration. From (4.11) we note that {a + l3)vo = avg + (3vf is constant, hence we dehne 
the equilibrium in such a way that the two displacement vectors i/g and i/f with respect to 
that equilibrium obey ayg + fdyj = 0 , and in particular we choose the equilibrium velocity 
to be vq = 0. We shall now dehne the relative displacement or dipole moment vector by 
i = yg- Vf-^ Hence (3.3) imply 


= (4.13a) 

= + (4.13b) 


where Pq is the common density of the two DM huids at equilibrium in the absence of external 
perturbations (far from any external mass), and P = Pq^ is the polarisation. The velocity 
helds of the two huids are 


/3 d| _ g d^ 

a + (3 dt ^ ^ a + (3 dt ^ 


(4.14) 


where d^/dt = dt$, -|- Vq ■ denotes the convective derivative, which reduces here to the 
ordinary derivative since Uo = 0. 

By inserting (4.13) into (4.5) we can solve for the polarisation resulting in 


P = -ApgWx'V.p. 


(4.15) 


®Here t/g and y/ denote the spatial components of the displacement four vectors and considered 

in Sec. 3, see (3.3)-(3.4). 
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Furthermore, by combining (4.14) and (4.11) and making use of the previous solution for the 
polarisation (4.15), we readily arrive at a simple harmonic oscillator describing plasma like 
oscillations aronnd equilibrium, namely 





(4.16) 




In this case the plasma frequency is given by 




a + (3 Tgpl 
(3 


(4.17) 


Next, consider the eqnation for the Newtonian potential Ug in the ordinary sector. Com¬ 
bining (4.10) with (4.4) we readily obtain 





(4.18) 


To recover the usual Newtonian limit we must impose, in geometrical units. 



(4.19) 


This condition being satished, the right-hand side of (4.18) can be rewritten with the help 
of (4.13) and the relation (4.12). We obtain an ordinary Poisson equation but modihed by 
polarisation effects. 


(4.20) 


V • VUg - Att— = . 

r„ 


We still have to show that the polarisation will be aligned with the gravitational held Vf/^, i.e. 
we grasp a mechanism of gravitational polarisation. This is a consequence of the constitutive 
relation (4.15) taken at the equilibrium point, neglecting plasma like oscillations. At this 
point we have Vf/g -|- = 0. Thus, 

P = ^Wx VUg. (4.21) 

Finally (4.20) with (4.21) takes exactly the form of the Bekenstein-Milgrom equation [64]. To 
recover the correct deep MOND regime we must impose that when A —)■ 0 

ivr/1 

1 _ 47r^ W;, = , (4.22) 

r] ao 

where Oq is the MOND acceleration scale. This is easily achieved with the choice 


>V(A) = A - |(a + + O [X^) , 

o 


(4.23) 
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together with hxing the constants Mefj and A to the values 

= A = |. (4.24) 

We thus recovered gravitational polarisation and the MOND equation (in agreement with 
the dielectric analogy of MOND [82]) as a natural consequence of bimetric gravity since it is 
made possible by the constraint equation (4.9) linking together the two metrics of bigravity. 
The polarisation mechanism is possible only if we can annihilate the gravitational force by 
some internal force, here chosen to be a vector held, and therefore assume a coupling between 
the two species of DM particles living in the g and / sectors. 

Let us recapitulate the various constraints we have found on the parameters in the original 
action (2.1). These are given by (4.12) for having a polarisation process, (4.19) for recovering 
the Poisson equation and (4.24) for having the correct deep MOND regime. In addition, we 
recall (3.8) which was imposed in order to be able to expand the two metrics around the same 
background. Strictly speaking the relation (3.8) is not necessary for the present calculation 
because we neglected the inhuence of the background, and (3.8) may be relaxed for some 
applications. Finally, we note that the result can be simplihed by absorbing the remaining 
constant Vg together with the sum a + f3 into the following redehnitions of the vector held 
and the effective metric: Ag —?• TgAg and —)■ (a + P)~‘^g'^- When this is done, and after 
redehning also the graviton mass m, we see that we can always choose Tg = 1 and a + (3 = 1 
without loss of generality. 

Finally the fully reduced form of the action (2.1) reads 


■^final \/ 'g 


[Ml 

— Rg— pbar ~ Pg 




1 - Ag,U^g 


+ \/—9eS 


dvr 87r ^ 



Pf 


1 + 




(4.25) 


where we have moved for convenience the mass currents = pgu'^ and Jj = ^Pfu'^ to the g 
and / sectors, where the effective metric g'^ of bigravity is given by (2.4) but with a + = 1, 

and where the kinetic term of the vector held is dehned by 


T’ 


irgiy tt 

fiiy 


W{A) = A- + O (T’2) . 
3 


(4.26a) 

(4.26b) 


In addition the coupling constants Mg and M| are constrained by (4.19). We can still further 
impose (3.8) if we insist that the two metrics g and / can be expanded around a common 

the previous model [87] we had to assume that some coupling constant e in the action tends to zero. 
With bimetric gravity and its nice potential interactions we see that no particular assumption is required. 
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background g = Finally the theory depends also on the graviton’s mass m, hopefully 
to be related to the observed cosmological constant A, and the MOND acceleration scale 
ag. Note that in such unihed approach between the cosmological constant and MOND, it is 
natural to expect that oq and A should have comparable orders of magnitudes, i.e. oq \/A 
which happens to be in very good agreement with observations. 


5 The decoupling limit 


In this section we would like to focus on the matter sector and address the question whether 
or not the interactions between the matter helds reintroduces the BD ghost [12], For this we 
derive the decoupling limit and pay special attention to the helicity-O mode. Before doing so, 
we hrst restore the broken gauge invariance by introducing the Stueckelberg helds in the / 
metric (with a,b = 0,1, 2, 3) 

ffj.h' —t ftMU = , (5.1) 


where the Stueckelberg helds can be further decomposed into their helicity-O vr and helicity-1 
A“ counterparts. 


= X 


mA°‘ f°‘°dh'x 


Ai 


Ai 


(5.2) 


As next we can take the decoupling limit by sending Mg, Mf —)■ oo, while keeping the scale 
Ag = MgW? constant. For our purpose, it is enough to keep track of the contributions of 
the helicity-O mode vr to the matter interactions. The BD ghost is hidden behind the higher 
derivative terms of the helicity-O interactions after using all of the covariant equations of 
motion and constraints. We will therefore neglect the contribution of the helicity-l held and 
assume and 


fIIV - t fgy - iVlMU Agy') 


(5.3) 


where stands for = dgdyii /In the remaining of this section all indices are raised 
and lowered with respect to the Minkowski metric g^y. The ehective metric in the decoupling 
limit corresponds to 

afy gfy = [(a + - f^Agy] ^ . (5.4) 

In [88] the required criteria for ghost freedom in the kinetic and potential terms were inves¬ 
tigated in detail and hence the new model was constructed in such a way that it fulhls these 
criteria. In the decoupling limit, the contribution to the equation of motion for the helicity-O 
held coming from the potential term is at most second order in derivative for the allowed 

simple choice satisfying all these requirement is a = /3 = t and Mg = M| = i.e. each coupling 
constant takes half the GR value. With this choice the plasma frequency (4.17) becomes 


(jj 


9 — 



in agreement with the finding of the previous model [87]. 
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potentials. Therefore, we can concentrate on the problematic terms coming from the matter 
interactions. Let ns for instance consider the conpling between and Afj,, thns^^ 


Llrna.t 


Pf \/ 9eS Afj,j j “I" ^ ^eff 1 


which we can also write as 


r — 

•^mat 


-\[A Pf{l + A^u‘t) + ■ 


(5.5) 


(5.6) 


Let ns hrst compnte the eqnation of motion with respect to the vector held, which is simply 
given by 



which can be also expressed as 




\/=/ . . 


7 — 


(5.8) 


As next we can compnte the contribntion of the matter to the eqnation of motion with respect 
to the helicity-O held, 

( (5£mat ^~9fa 


^-^mat ^fpa j ^ '-'fi'-'ip 

StT * ^ r-i-r 


Ai V 57 


Ai \ 


d/,du 

Ai 

where we made nse of 

rjnflU _ 



V7 

Applying one of the derivatives we have 
^^mat dy f _ 


j rriUU 2 (5/lniat 

and 1 ^ = , ——^ . 

\/-ges ^gfu 


p p> 




5ti ■' f 

+ / 3 a„( v^ri'')[(a + , 3 ) 5 “ - , 3 n“] - 9 „n 

The eqnation of motion for the vector held (5.8) can also be written as 

'HP TSeff — T 


(5.10) 


(5.11) 


(5.12) 


^For simplicity we will drop the constants A, and M^g and assume that yV{X) = X in the following. 
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which we can use to express 


K-’^Z = +r;r r- = jpsj-.,. (5.13) 

hence we have 

a„ {v^«ZZ) = sS r., - \/^r^-*rr. (5,i4) 

Furthermore, the conservation law [see (2.15)] gives 

V^^K-TZ + = 0 ■ (5.16) 

Thus, we have further 

d, U^fTz'\ = - '/Qr^Tf. (5.16) 

Using the equation of motion for the vector held and the conservation equation, the equation 
for the helicity-O becomes 

^ = -4 { - 'fZrZTf'j (s; - no - ( 5 . 17 ) 

+ P - V^rj“r£) [(a + i3)s; - . 

which we can rewrite as 

- n;;) - 'fZTpn-^ 

+ aZ^dp/Za^ZZ + m - /sna - ZZzZtzP''^] , ( 5 . 18 ) 

where we have introduced 

KO = rV(5;-no + a"n„,, ( 5 . 19 a) 

iiZ = rjl" [(“ + ZK -K\+ . (6. i9b) 

Using the fact that the Christoffel symbol with respect to the / metric is given by 

rd = -'r(K - n09„n„,, ( 5 . 20 ) 

and further taking into account the following relation, 

(5.21) 
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we see immediately that = 0. Similarly, using the dehnition of the Christoffel symbol 
with respect to the ^eff metric, 

r;r = -&(» + m - , (5.22) 

and the relation 

[{a + l3)S;-mf\S[(a + m-f>K]=^r, (5-23) 

we can show that also = 0. Thus, the contribution of the matter part to the equation of 
motion for the helicity-0 mode simplihes to 

^ - np++^))y - ,3np}, 

(5.24) 

Using (5.21) this can be equally written as 

^ - n"")-' + [(a + /3)ir - /sn-"]I , (5.25) 

where = \/—gesjj is the coordinate current and represents the matter degrees of freedom 
independent of the metric, while is also independent of the metric. We can now also apply 
the remaining derivative in front, 

- np + 

+pdAX&A [(o + P)s; - pn;] - . (5.26) 

The contribution of the matter interactions to the equation of motion for the tt held contains 
higher derivative terms which we are not able to remove by invoking the covariant equations 
of motion for the vector held and the dark matter particle. Hence, this rehects the presence 
of a ghostly degree of freedom in the decoupling limit through the matter interaction term. 
The scale A 3 is not the cut-oh scale of this theory anymore, but rather given by the scale of 
the mass of the introduced ghost 772 , 30 - For the precise plasma-like background solution that 
we considered here, the contributions of the matter interactions and the internal vector held 
will be at least third order in perturbation at the level of the action. Hence the ghost would 
enter at non-linear order in perturbations, and does not show up in our linear perturbation 
analysis of the gravitational polarisation in Sec. 4. This indicates that the mass of the ghost 
772 bd could be large. 

6 Conclusions 

In this work we followed the philosophy of combining the diherent approaches used to tackle 
the cosmological constant and dark matter problems. The origin of this dark sector consti¬ 
tutes one of the most challenging puzzle of contemporary physics. Here we want essentially to 


6C 


mat 
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consider them on the same footing. The standard model of cosmology A-CDM, despites many 
observational snccesses, fails to explain the observed tiny value of the cosmological constant 
in the presence of large quantum corrections using standard quantum held theory techniques. 
This has initiated the development of IR modihcations of GR, like massive gravity and bigrav¬ 
ity. On the other hand, the model A-CDM does not account for many observations of dark 
matter at galactic scales, being unable to explain without hne-tuning the tight correlations be¬ 
tween the dark and luminous matter in galaxy halos. Rather, the dark matter phenomenology 
at galactic scales is in good agreement with MOND [61-63]. 

In the present work we aimed at addressing these two motivations under the same um¬ 
brella using a common framework, namely the one of bigravity. An important clue in this 
respect is the fact that the MOND acceleration scale is of the order of the cosmological con¬ 
stant, Oq ~ \/A. An additional hope was to be able to promote the MOND formula into a 
decent relativistic theory in the context of bigravity. For this purpose we followed tightly the 
same ingredients as in the model proposed in [88], with two species of dark matter particles 
coupled to the two metrics of bigravity respectively, and linked via an internal vector held. 

This paper was dedicated to explore the theoretical and phenomenological consistency 
of this model and verify that it is capable to recover the MOND phenomenology on galactic 
scales. We hrst worked out the covariant held equations with a special emphasis on the 
contributions of the vector held and the dark matter particles. Because of the interaction 
term between the dark matter particles and the vector held, the stress-energy tensors are not 
conserved separately, but rather a combination of them, giving rise to a global conservation 
law. The divergence of the stress-energy tensor with respect to the ehective metric Qes to which 
is coupled the vector held [15, 17], is given by the interaction of the dark matter particles with 
the vector held, which has important consequences for the polarisation mechanism but also 
for our decoupling limit analysis. 

As next we computed the linear held equations around a de Sitter background where the 
mass term shall play the role of the cosmological constant on large scales. On small scales 
where the post-Newtonian limit applies, the de Sitter background can be approximated by a 
hat Minkowski background. Considering a small perturbation of the Minkowski metric and 
computing the Newtonian limit, we were able to show that the polarisation mechanism works 
successfully and recovers the MOND phenomenology on galactic scales. We hnd that this 
polarisation mechanism is a natural consequence of bimetric gravity since it is made possible 
by the constraint equation (linearized Bianchi identity) linking together the two metrics of 
bigravity. In addition it strongly relies on the internal vector held generated by the coupling 
between the two species of DM particles living in the two metrics, as it is alive due to an 
appropriate cancelation between the gravitational force and the internal vector force. 

We then studied the decoupling limit of the matter sector in details. We found that 
precisely because of this coupling to the vector held (whose importance is probably priceless 
for the polarisation to work) a ghost is reintroduced in the dark matter sector. For the 
particular solution we have for the dark matter the ghost is higher order in the perturbative 
expansion. The exact mass of the ghost shall be studied in a future work. 
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